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Abstract. The effects of an aperiodic order or a random disorder on phase transitions 
in statistical mechanics are discussed. A heuristic relevance criterion based on scaling 
arguments as well as specific results for Ising models with random disorder or certain 
kinds of aperiodic order are reviewed. In particular, this includes an exact real-space 
renormalization treatment of the Ising quantum chains with coupling constants modu- 
lated according to substitution sequences, related to a two-dimensional classical Ising 
model with layered disorder. 



1 Introduction 

Equilibrium statistical mechanics bridges the gap between the fundamental laws 
of physics at microscopic scales on the one hand and the macroscopic phenom- 
ena that we may experience in every-day life on the other hand jy. Its main 
achievement is the understanding of the behaviour of macroscopic systems on 
the basis of the microscopic constituents and their mutual interactions. 

One of the most fascinating topics in this field concerns phase transitions, 
i.e., a qualitative change observed in some physical properties of a macroscopic 
system after a slight variation of certain parameters as, for instance, the tem- 
perature or an external magnetic field. Here, statistical mechanics reveals its 
full power by being able to explain how even short-range interactions at the mi- 
croscopic scale may lead to such violent cooperative behaviour. In mathematical 
terms, a phase transition corresponds to a non-analyticity of the thermodynamic 
potential that describes the equilibrium state of the system as a function of the 
external parameters. In essence, one may regard it as a battle between energetic 
and entropic terms, the first driving the system towards an ordered phase at low 
temperature, the second favouring a disordered state at higher temperatures. 
However, this is a somewhat simplistic view of a very complex phenomenon that 
should not be overstressed. As an example that teaches one to be cautious one 
might think of the entropy-driven phase transitions observed in non-interacting 
hard-core systems || that, counter-intuitively, may lead to an ordered state that 
has a higher entropy than the disordered state. 

While the general picture is rather well understood, it turns out to be much 
more difficult to calculate the properties of a specific macroscopic system from its 
microscopic interactions explicitly. In fact, only for a few simple models complete 
solutions are known |3| . Among these is the paradigm of a ferromagnet, the 
celebrated Ising model [ |,|) , after Onsager || accomplished to calculate the free 



2 



Uwe Grimm 



energy for the two-dimensional Ising model without an external magnetic field. 
However, also this extremely simplistic model has, so far, defied any attempts 
at an analytic solution in higher dimensions, or even in two dimensions in the 
presence of an external magnetic field. Recent progress from complexity theory 
might even support the widespread belief of its general unsolvability. 

Moreover, the models that have been solved analytically (or 'exactly', which 
usually does not mean rigorously in the mathematical sense) are almost exclu- 
sively tied to regular periodic structures, as for instance the Ising model on a 
square lattice or some other regular planar lattices. However, Nature is never 
completely regular, and this immediately spurs the question what happens if 
the regular periodic structure is disturbed by some randomness, mimicking a 
more realistic crystalline substance, or if it even is given up completely, resulting 
in either aperiodically ordered, or disordered structures. While the effects of a 
random disorder have been studied for many decades, aperiodically ordered sys- 
tems entered the scene after the experimental discovery of aperiodically ordered 
solids known as quasicrystals in the 1980s, see for a recent compilation of 
introductory lectures on this issue. 

The situation is, however, not as hopeless as it may seem on first view. One 
of the things that come at our rescue is the important concept of universality 
of critical phenomena (|]. According to the universality hypothesis, the critical 
behaviour of a system does not depend on details of the microscopic interactions 
of the model, but merely on a number of general properties such as the dimension 
of space and the symmetries of the model, see also the contributions P Jl0|[ll] , |l2| 
in this volume. This means that second-order or 'smooth' phase transitions, 
which are characterized by an infinite correlation length, can be classified into 
universality classes of systems showing the same physical behaviour at and close 
to the critical point. In particular, this implies that the critical exponents that 
describe the singularities of physical quantities upon approaching the critical 
point only depend on the universality class of a model. The origin of universality 
and the scaling relations between different critical exponents can be understood 
in terms of the so-called renormalization group []l3|]9f| because critical points 
correspond to fixed points under the renormalization flow. Of course, in reality 
things are not quite as simple, and the situation may, in general, turn out to be 
more complex. Therefore, such arguments should be taken with a grain of salt. 

In two dimensions, we have another powerful tool at our disposal. There ex- 
ists an intimate relation between two-dimensional statistical models at criticality 
and conformal field theory ]l4| that explains the surprising fact that critical ex- 
ponents in many cases are found to be rational numbers, and furthermore allows 
the computation of correlation functions. This yields an alternative, and more 
substantiated, explanation of the universality classes and a partial classification 
of the possible types of critical behaviour in two dimensions. Such strong, albeit 
by no means complete, results are limited to the two-dimensional case due to 
the particular properties of the two-dimensional conformal symmetry group. 

Nevertheless, because of universality, one may expect a certain robustness 
of critical point properties, i.e., introducing disorder or aperiodic order may not 
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affect the critical behaviour of some systems at all. The obvious question one 
would like to answer is, given an aperiodic or disordered model, does it show the 
same critical behaviour as its perfectly ordered fellow or not? And, if not, how 
does it change? 

In what follows, we shall concentrate on the example of the Ising model, 
simply because it is the most frequently studied and, besides numerous numerical 
treatments, also permits analytic solutions. After a brief introduction to the 
Ising model and its critical properties in the subsequent section, we discuss 
a rather general, albeit heuristic, relevance criterion that is based on scaling 
arguments in section ||. Thereafter, in section ^, we give a short overview on 
the various types of disordered Ising models considered in the literature and the 
main results concerning their critical behaviour. Next we turn our attention to 
aperiodically ordered Ising models which are the topic of section [|. Finally, the 
article concludes with a brief summary in section [| 



2 The Ising Model 

There have been numerous review articles and even entire books jl5| devoted 
to the Ising model, and many of these give en extensive introduction into the 
history of this simple, yet fascinating model. Also most textbooks on statistical 
mechanics treat at least the one-dimensional Ising model. For a mathematically 



rigorous approach, see, e.g., |16 . Here, we shall restrict ourselves to the most 
important facts, the interested reader is referred to the literature for a more com- 
plete account of the historical background, for instance in |17|. Also the recent 



review |18 on aperiodic Ising models contains a brief historical intreoduction, 



and an extensive list of references. 



2.1 The Square-Lattice Ising Model 

Originally, the Ising model was devised by Lenz as a simplistic toy model of 
a ferromagnet The one-dimensional model was later solved by Ising a 
student of Lenz, who obtained the disappointing result that the model does not 
show a phase transition at a non-zero temperature. He supposed that this should 
hold in higher dimensions as well, but, fortunately, it did not take too long until 
he was proven wrong - otherwise the model carrying his name would not have 
attracted the interest of physicists for almost a century. 

So, how is the model defined? The local magnetic moments of the solid are 
modeled by 'Ising spins' on the sites of a regular lattice, which can only take two 
values, ±1. The ferromagnetic interaction between these spins is given by a sim- 
ple coupling of neighbours on the lattice. Let us concentrate on the square-lattice 
case for definiteness. A configuration a — {p-^ \i o"i 2i • • • i &n m 

) e {±1} NM of 

spins <r ■ k on a N x M rectangular section of the lattice is assigned an energy 

N M 

E (°) XX ( K(T J.k°3 + l,k + L °j,k°j,k+1 + H07 jk ) (1) 

3 = \ fc=l 
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Fig. 1. The square-lattice Ising model 



where we may assume certain boundary conditions such as, for instance, periodic 

Ka/+i = a N+i.k = or frce ( CT j,M+i = °jv+i,k = °) boun dary condi- 
tions. Here, K > and L > denote the ferromagnetic coupling constants in 
horizontal and vertical direction, respectively, see Fig. |l|, and H is the magnetic 
field. 

The partition function, a sum of Boltzmann factors over all configurations a, 
is given by 

Z N , M (T,H) = ]Texp[-£:(a)/fc B T] (2) 

a 

where T denotes the temperature and fc B is Boltzmann's constant. It can also 
be expressed in terms of a so-called transfer matrix 

Z NM (T, H) = Tr[U N (T, H) m ] (3) 

where U N is a 2 N x 2 N matrix whose elements are Boltzmann weights of the 
configuration corresponding to a single row of the lattice, and the trace is due 
to the periodic boundary conditions that we assumed. Thus, the trace performs 
the sum over the configuration (a 1 1 , <t 2 1( . . . , a N x ) of the first row, whereas the 
summation on the configurations {a 1 fe , a 2 ^, • • • , cr N fc ), 2 < k < N, along the 
other rows are hidden in the matrix product. 
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2.2 Phase Transitions and Critical Exponents 

In the thermodynamic limit of an infinite system, the free energy per site 

is an analytic function of the parameters T and H, except at the phase transi- 
tion. The transition is of first order if already a first derivative of the free energy 
is discontinuous. Otherwise, the transition is termed a second- or higher-order 
transition. Examples of first-order transitions are melting of ice or the evapo- 
ration of water, where one has to supply energy, the latent heat, to the system 
without changing the temperature. The Curie point of a ferromagnet, or the 
critical point of water at the end of the first-order transition line between liquid 
and gas, are examples of critical points. 

Correspondingly, one encounters singularities in thermodynamic quantities 
that can be obtained as derivatives of the free energy such as the specific heat 

at constant magnetic field H, the spontaneous magnetization 

m (r)=lkm(T,F), m(T, H) = 9f f^ H) , (6) 
or the magnetic susceptibility 

Xt{H) ~ dH 8IP- (7) 

For a second-order phase transition, on which we shall now concentrate, these 
singularities are characterized by so-called critical exponents describing power- 
law singularities on approaching the critical point. For magnetic systems such 
as the Ising model, some zero-field (H = 0) critical exponents are 

c (T) (8) 

m o( T ) ~ (-if (for t < 0, m (T) = for t > 0) (9) 

Xr (0) ~ \tr (10) 

as t — > 0, where t = (T — T c )/T c is a 'reduced temperature' parametrizing the 
distance to the critical temperature T = T c , see also Q. The power- laws are 
related to a particular behaviour of correlation G(r, r') at the critical point. In- 
stead of a usual exponential decay with the distance, G(r, r') ~ exp(— \r — r'|/£), 
with a characteristic correlation length £, correlation functions at a critical point 
show a power-law decay G(r,r') ~ \r — r'\~ x ). In other words, the correlation 
length £ diverges as 

e~i*r (ii) 

as the critical point is approached. Of course, several other exponents can be 
considered, and there also exist phase transitions where the exponents differ 
when the critical point is approached from the ordered phase t < or from the 
disordered phase t > 0. 
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2.3 The Ising Model Phase Transition 

But as mentioned above, it was not clear from the beginning that the Ising 
model has a phase transition at all. However, already in 1936, Peierls JjiJ gave 
an argument that showed the existence of a phase transition in two or more 
dimensions. The argument is simple, and indeed also holds for rather general 
models, including Ising models with ferromagnetic random couplings. 

Essentially, the argument goes as follows. At T = 0, the system is in one of 
its magnetized ground states, where all spins are either +1 or —1. The states 
that the system can explore at a small temperature differ from the ground state 
by small islands of turned spins, and their energy is given by the surface of these 
islands where neighbouring spins differ. So, even at small finite temperatures, 
the system is magnetized, and there must be a transition to a disordered high- 
temperature state at a non-zero temperature. Indeed, this argument breaks done 
in one dimension, because a flip of the local spins along an entire infinite half- line 
affects a single bond only, thus even an arbitrarily small temperature suffices to 
destroy the magnetic order, see |2(J for a more precise discussion. 

The critical temperature for the zero-field Ising model was first determined 
by Kramers and Wannier [^l] in 1941 under the assumption that there is only 
a single phase transition. In this case, the transition has to take place at a fixed 
point of a duality transformation that maps low- and high-temperature phases 
onto each other. For the square-lattice Ising model, the critical temperature T c 
is given by 

smh(2K/k B T c ) sinh{2L/k B T c ) = 1 , H = , (12) 

Only three years later, Onsager's celebrated solution of the zero-field square- 
lattice Ising model || revealed a logarithmic singularity of the specific heat, 

co(t) ~ ln|l/t| , (13) 

i.e., the corresponding critical exponent a — 0. The correlation length diverges at 
the critical point with an exponent v = 1 . The magnetization exponent f3 = 1/8 
was later calculated by Yang |^2| by means of a perturbative approach. In Fig. ||, 
the specific heat and the spontaneous magnetization for the zero-field Ising model 
are shown. 




Fig. 2. The specific heat Cq(T) and the spontaneous magnetization mo(T) of the 
square-lattice Ising model 
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2.4 The Ising Quantum Chain 

Before we can, finally, move on to discuss disordered models, we need to intro- 
duce a one-dimensional quantum version of the Ising model |23|| that is closely 
related to the two-dimensional classical Ising model, see also |12| . Consider the 
Hamiltonian 

/ n n \ 

H = — 2 (E°W+i +A E°i) ( 14 ) 

on the Hilbert space 0^ =1 C 2 S C 2 . Here, the local spin operators are defined 



as 

X, z 



(g)I(g)...(g)I(g)I (15) 



j — X factors N—j factors 

with the standard Pauli matrices 

1 oj ' ° = 0-1 J ' I= lo 1 



(16) 



Frequently, a different basis is used, where the coupling term involves the prod- 
ucts and the 'transversal field' term has the form J2f=i a j Pjl - For the 
homogeneous Hamiltonian ( |i~4| ) , the spectrum is known completely J24| . 

The Hamiltonian ( [Ti] ) can be obtained as an anisotropic limit of the transfer 
matrix U N (||) of the zero-field (H — 0) square-lattice Ising model (Q) , where the 
coupling along the chain direction K — > and the coupling in the 'time' direction 
L — ► oo, while 0K exp(2/3L) is kept fixed p5j ]. Consequently, the 'transversal 
field' parameter A in the quantum chain acts analogously to the temperature in 
the classical model; and the Ising quantum chain at zero temperature displays 
a quantum phase transition Jl2[ , i.e., a change in the ground-state properties, 
at the critical value A c = 1 which belongs to the same universality class as 
the critical point of the classical square-lattice Ising model. In particular, the 
ground-state energy per site, —E /N, of the Hamiltonian ( p^ ) corresponds to 
the free energy of the classical Ising model, and the energy gap E l — E Q between 
the ground state and the first excited state corresponds to the inverse t;^ 1 of 
the correlation length along the chain, i.e., the energy gap vanishes at criticality. 
According to finite-size scaling |^7|, the gap vanishes as E 1 — E ~ N~ z , where 
z = denotes the ratio of the correlation length exponents along the chain 

and in the 'time' direction. In spite of the anisotropic limit, at criticality the 
homogeneous system (|l4| ) behaves isotropically, with ^ = i/ ± = 1 as in the 
classical square-lattice Ising model. However, as we shall see below, this may 
be different if one considers disordered or aperiodically ordered Ising quantum 
chains, which then correspond to classical two-dimensional Ising models with a 
layered disorder or aperiodicity. 



2.5 Effects of Disorder: Some General Remarks 

Phase transitions are cooperative phenomena leading to singularities in physical 
quantities such as the specific heat or the magnetic susceptibility. From a naive 
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point of view, any disorder in a system, if it affects the phase transition at 
all, will tend to weaken critical singularities - it is hard to imagine how it could 
possibly do the opposite. Thus, one might expect that a first-order transition in a 
perfectly ordered system may be weakened to a higher-order phase transition in a 
disordered or an aperiodically ordered system, and indeed there exist examples 
where such behaviour has been observed |26[ |. For a second-order transition, 
disorder may change the critical exponents, or the singularities may be weakened 
to higher-order singularities, or even washed out completely so that the phase 
transition disappears. 

So, given a specific model, how can one estimate the effect of a certain type 
of disorder on the system? A quite simple answer to this question is provided 
by a heuristic relevance criterion based on scaling arguments. Although these 
arguments rely on several assumptions that may not always be fulfilled, they 
have proven to be rather successful in predicting the correct critical behaviour. 



3 Heuristic Scaling Arguments 

The relevance criterion that we are going to discuss now had first been put for- 
ward be Harris in 1974 |2^] for disordered Ising models. Later, it was generalized 
by Luck J2l| to the case of aperiodically ordered Ising models. The argumenta- 
tion presented here closely follows the discussion of the Harris-Luck criterion in 
[ p7| , where one can also find further applications of the criterion. 

We restrict ourselves to models of Ising type with purely ferromagnetic cou- 
plings; the case of frustration in disordered systems or spin glasses, see |$0|, is 
out of reach and thus not discussed here. Furthermore, a certain homogeneity 
of the distribution of the ferromagnetic coupling constants £ • fc > is needed, 
such that, for instance, the mean coupling e = e is well defined, which we 
shall consider as the coupling of our unperturbed reference system denoted by 
a lower index 0. We characterize the distribution of coupling constants by their 
fluctuations, expressed in terms of the deviations from the mean coupling. More 
precisely, consider an approximately spherical volume V located somewhere in 
our infinite system, see Fig. |3[ Within V , the mean coupling is 

£ ~ v= ^k,^ e ** (17) 

where N v denotes the number of spins in V, and n is the mean coordination 
number, i.e., the averaged number of neighbours in the system. The deviation 
of the accumulated couplings from the mean is given by 

As a measure of the fluctuations, we consider the asymptotic behaviour of the 
standard deviation A G for large volumes V, 

A Gv ~ (AV)" , (19) 
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Fig. 3. A volume V within the system and the associated local quantities 

thus defining the fluctuation exponent uj which we assume to be well defined. 
Note that we did not specify the type of disorder in the system, we may consider 
both models with varying coupling constants on a regular lattice, or models that 
are defined on topologically disordered graphs, for instance random graphs or 
periodic lattices with randomly distributed vacancies. Another example we shall 
encounter below concerns Ising models that are defined on aperiodic graphs, 
which also yield fluctuations albeit there is no randomness involved. 

Now comes the crucial intuitive step. Clearly, the critical temperature T c 
depends on the average coupling constant e = e Q . Thus, due to the fluctuations 
in the couplings e - fc , the system may locally, i.e., within a volume V, correspond 
to a different 'local critical temperature' T c y that, for sufficiently small disorder, 
may be expected to depend linearly on the local average coupling T c v ~ e v . 
One may think of Tj as the critical temperature of a system that essentially 
looks the same in any ball of volume V. Correspondingly, we may define a 'local 
reduced temperature' t v . 

If the disorder is irrelevant for the critical behaviour, the shift in the local 
reduced temperature &t = t v — t has to vanish at criticality, i.e., St — > as 
t — > 0. The relevance criterion now follows from this consistency requirement 
by using scaling relations involving the critical exponent which in the case of 
irrelevant disorder should coincide with the exponent v of the pure system. The 
argument proceeds as follows. The volume V of spins that are correlated is of 
the order V ~ £ d , where £ denotes the correlation length, and d the dimension 
of the system. Then, from ( jl9| ) and our assumption T c y <~ e v , we deduce 



where v is the correlation exponent of the disordered system. In the irrelevant 
scenario v = v , this yields 




(20) 



bt/t ~ \t\ 



with <P = 1 - dv Q {l - to) . 



(21) 
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Consequently, the disorder is irrelevant if the crossover exponent < 0. 

We can summarize the Harris-Luck relevance criterion as follows p8|j29|p7| |. 
A modulation in d m < d space dimensions of a ferromagnetic <i-dimensional Ising 
model with correlation exponent v is relevant for lo > lo c , marginal for lo = lo c , 
and irrelevant for lo < lo c1 where 

«„ = 1-^ (22) 

and lo denotes the fluctuation exponent (|l9|). 

For the two-dimensional Ising model, we have v = 1. Thus, for a model 
with layered disorder, i.e., d m = 1, we obtain lo c = 0. In this case, any divergent 
fluctuation is relevant. For planar disorder, i.e., d m = 2, the marginal value 
is lo c = 1/2. Thus, for randomly distributed couplings, which correspond to 
lo = 1/2, one is precisely in the marginal situation where the criterion does 
not give a definite prediction. Generally, one might expect to find logarithmic 
corrections to scaling in the marginal case, and this is indeed the case in the 
two-dimensional random-bond Ising model |34|j35| , |3^ , ^7| . 

In the original formulation by Harris, the criterion was expressed in terms 
of the specific heat exponent a in place of the correlation exponent v. The two 
exponents are related by the hyperscaling relation dv = 2 — a 0, which holds 
for the two-dimensional Ising model, which again corresponds to the marginal 
case a = in the case of random disorder [E8|. 



4 Random Ising Models 

Traditionally, randomness was introduced into Ising models to describe magnetic 
behaviour of alloys where some magnetic atoms were replaced by non-magnetic 
atoms. Therefore, such models are termed dilute Ising models, and, in the sim- 
plest case, just correspond to a regular Ising model where some spins, or some 
bonds, have been removed, see the reviews ||l],[52|,[33| . Clearly, there is some re- 
lation to percolation ]Ic|| , because one at least needs infinite clusters in order 
to support a phase transition. Somewhat more general are random-bond Ising 
models, where the coupling constants are chosen randomly, often from a bimodal 
distribution. 

In these systems, the disorder is regarded as 'frozen' or 'quenched', which 
means that it is static. Therefore, one has to average on the free energy on 
the possible realizations, which is extremely difficult |3(J. Nevertheless, some 
analytical results have been obtained. As shown above, we expect the two- 
dimensional system to be marginal with respect to random disorder. Using the 
replica method, Dotsenko and Dotsenko (34| found a double logarithmic singu- 
larity in the specific heat 

c(t)~]n[l+0]n(l/|t|)], (23) 
in place of the pure Ising behaviour (|l^) , and 
m(t) ~ cxp{-a[lnln(l/K|)] 2 /2} , *(*) ~ t~ 2 exp{-a[lnln(l/|i|)] 2 } , (24) 
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corresponding to critical exponents (3 = and 7 = 2, which differ from the pure 
Ising values. However, by means of bosonization techniques and conformal field 
theory, Shalaev, Shankar, and Ludwig 35 3(]j37j arrived at different results 



m(t) ~ t a -/ 8 0n|*l) _x/1 ° , X(t) -^ 7/4 (ln|t|) 7/8 , (25) 

for the magnetization and the susceptibility, while the result for the specific 
heat coincides with (j^TJ). Here, the critical exponents keep their pure Ising val- 
ues (3 = 1/8 and 7 = 7/4, but logarithmic corrections show up. Due to the slow 
variation of the logarithm, it turns out to be very difficult to obtain reliable 
results from numerical calculations, because the logarithmic dependence on the 
reduced temperature translated to a logarithmic dependence on the system size, 
and it is hopeless to distinguish a simple logarithmic from a double logarith- 
mic from finite-size calculations. Still, the majority of the numerous finite-size 
scaling studies based on Monte Carlo simulations or transfer matrix calculations 
favour the latter result ( p5| ) , see 1 31 3^,|33| and references therein, which was also 



substantiated by a recent series expansion investigation 

The correlation exponent of the three-dimensional pure Ising model v w 0.63, 
thus a > 0, and disorder should be relevant in that case. However, a recent result 
indicates that for purely topological disorder introduced by considering an Ising 
model on a three-dimensional random graph the critical behaviour stays the 
same as for the regular cubic lattice ]39| ]. 

In recent years, there has been an increasing activity in the mathematical 
literature on stochastic Ising model, see, e.g., ]4fj|j4l| , ^2t , leading to a number of 
rather general rigorous results. For instance, for the random-field Ising model, 
it can be shown rigorously that for dimensions d < 2 there exists a unique 
Gibbs state, thus there is no phase transition in the two-dimensional random- 
field model, whereas there is a phase transition in d > 3. For ferromagnetic 
random-bond Ising models on certain periodic planar graphs, it was proven that 
at most two extremal Gibbs states exist, so there are no more than two phases, 
see [E3| and references therein. 



5 Aperiodic Ising Models 

In the remainder of this article, we shall concentrate on the case of aperiod- 
ically ordered Ising modelsj TJ] . To this end, we need to introduce aperiodic 



sequences and tilings, see |44] for a recent review of the mathematical aspects. 
In the theoretical description of quasicrystalline structures 0, aperiodic tilings 
constitute the the analogues of periodic lattices in conventional crystallogra- 
phy. For the sake of space, we cannot go into much detail here. The interested 
reader is referred to the Mathematica |45| program packages accompanying this 
article, which give an insight into the construction and properties of one- and 
two-dimensional aperiodic tilings. The different methods employed and the cor- 
responding Mathematica routines were described in detail in [49]. 
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5.1 Aperiodic Tilings on the Computer 

Aperiodic sequences are commonly constructed by means of substitution rules. 
As examples, consider an alphabet of two letters a and b and rules 

C «: (* = 1,2,...) (26) 

that replace the single letters a and b by words Wa = ab and io^ = a k , 
the word consisting of k letters a, in the two-letter alphabet {a, b}. We restrict 
ourselves to primitive substitutions, i.e., after a finite number of iterations the 
words obtained from the basic letters should contain all letters. Corresponding 
semi-infinite substitution sequences w^) are then obtained by as fixed points of 
the substitution rules by an iterated application = g i - k \w^}_ 1 ) of the rules 

( p6|) on some initial word, say Wq = a. The substitution rules (|2^) yield some 
prominent examples of aperiodic sequences, for instance the Fibonacci sequence 
abaababa . . . for k = 1, the so-called periodic-doubling sequence abaaabab ... for 
k = 2, and a sequence known as the binary non-Pisot sequence abaaaababab . . . 
for k = 3. 

Many properties of the substitution sequences can conveniently be calculated 
form the associated substitution matrices 



M» = (;*), ^>_l±^H, (27) 
whose elements just count the number of letters a and b in their substitutes 



and , respectively. The largest eigenvalue of determines the 

asymptotic growth of the sequence in one iteration step, while the entries of 
the corresponding eigenvector, suitably normalized, determine the frequencies 
p a = 1 — p b = 1 — 2/(3 + V4fc + 1) of the letters a and b in . The second 
eigenvalue measures the fluctuations of the letter frequencies 

gW(N) = # a [w£\ N ] - Pa N , h^(N) = max \g^(M)\ ~ N" (k> , (28) 

where # a [wco \n] denotes the number of letters a in the first N letters of the 

limit word w£, ■ Thus, g( k \N) measures te deviation from the mean p a N, and 
we use the maximal deviation h^(N) to define the fluctuation exponent u>^ k ' 
which enters the Harris-Luck criterion, compare ([l9|). It is given by 

(fc) = H^| W (D = _ 1 W C3) =0 , c( 3 )«0.317>0; (29) 
InAf 

and thus, according to the Harris-Luck criterion, a layered aperiodicity intro- 
duced into the Ising model according to these substitution rules is irrelevant for 
k = 1, marginal for k = 2, and relevant for k — 3. So, contrary to what one might 
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Fig. 4. Eightfold patch of the Ammann-Beenker tiling 



have expected, deterministic substitution sequences, albeit being non-random, 
actually provide examples for all interesting classes of 'disorder'. 

The Mathematica program FibonacciChain .m that is included in the pro- 
gram package offers the reader the opportunity to gain experience with substitu- 
tion sequences and their properties. In addition, it also shows how quasiperiodic 
examples, such as the Fibonacci chain, can be derived from a projection of the 
two-dimensional square lattice. This approach, known as the cut-and-project 
method, can be generalized to higher dimensions to obtain quasiperiodic tilings 
(or in mathematical terminology, model sets [fti] ) with interesting symmetry 
properties. Roughly, a certain part of the high-dimensional periodic lattice is 
projected onto a suitably chosen subspace. As an example, Fig. || shows a cen- 
tral patch of an eightfold symmetric planar tiling obtained from projection of the 
hypercubic lattice Z 4 . The tiling consists of squares and rhombi and is known 
as the octagonal or the Ammann-Beenker tiling |^7j . In addition, this geometric 
structure also admits an inflation/deflation symmetry that is analogous to the 
substitution rule for the one-dimensional sequences. It consists of a dissection 
of the two basic tiles into copies of itself, such that the resulting tiling is, apart 
from an overall scaling, invariant under the procedure. 

Both approaches are used in the Mathematica program OctagonalTiling.m 
| fl6| . The remaining programs ChairTiling .m and SphinxTiling.m deal with 
two further planar examples of inflation tilings, while GridMethod .m introduces 
a variant of the cut-and-project scheme Efl. The package PenrosePuzzle .m 
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employs yet another method to construct the most famous among the zoo of 
quasiperiodic tilings, Penrose's tiling of the plane. Here, the two basic rhombi 
are marked by certain arrow decorations on their edges, and it is the task to fill 
the entire plane with tiles without holes or overlaps, and without violating the 
matching rules that require that arrow decorations on adjacent edges have to 
match. The reader should be warned that this, albeit seemingly simplest, puzzle 
method is not suited to construct a large patch of the tiling because the matching 
rules do not uniquely specify how to add tiles, and, in particular, do not prevent 
implicit mistakes which will only be revealed when, at a certain stage, no further 
legal additions of tiles are possible. But, if you indeed managed to fill the entire 
plans without violating the rules, then the resulting tiling is a Penrose tiling 
with all the other 'magic' properties such as an inflation/deflation symmetry 

5.2 Ising Models on Planar Aperiodic Graphs 

In what follows, we shall consider planar Ising models defined on quasiperiodic 
graphs, and, subsequently, Ising quantum chains with coupling constants varying 
according to substitution sequences. The latter correspond to two-dimensional 
classical Ising models with a layered aperiodicity, thus d m — 1 in Harris-Luck 
criterion, and lu c = 0. 

It may be surprising, but it is indeed possible to construct planar aperiodic 
Ising models that are exactly solvable ^8|,[l8| in the sense of commuting transfer 
matrices ||. For such systems, the coupling constants have to be chosen in a 
particular way to ensure the integrability |4S[] . As a consequence, they have the 
rather unusual property that the free energy does not depend on the actual 
distribution of coupling constant, but only on the frequencies of the various 
couplings. Another evidence of the particular arrangement is the fact that the 
local magnetization, i.e., the expectation value of a spin at a certain site, does not 
depend on the site, but is uniform on the entire system. Therefore, the geometric 
arrangement of couplings does not play a role, and the critical behaviour is 
always the same as in the pure system. Nevertheless, these models are interesting 
because they provide counter-examples to the Harris-Luck criterion, because the 
previous statement also holds for arrangements with strong fluctuations that 
should have been relevant. However, the solvable examples are certainly very 
special, and possess some hidden underlying symmetry; and one should expect 
that the criterion holds for any generic distribution of coupling constants. Still, 
it reveals the limited predictive power of such criteria. 

Ising models on quasiperiodic graphs, in particular the Penrose tiling, have 
been thoroughly investigated by means of Monte Carlo simulations and by ap- 
proximative renormalization group treatments, see [lq] and the literature cited 
therein. Usually, the coupling is taken uniform along the bonds of the graph, 
so the fluctuations arise solely from the locally differing coordination numbers. 
This type of planar aperiodicity is irrelevant according to the Harris-Luck crite- 
rion, because the fluctuations for quasiperiodic cut-and-project sets, such as the 
Penrose tiling or the Ammann-Beenker tiling shown in Fig. |J, are small. Essen- 
tially, this is due to the fact that the quasiperiodic tiling is the projection of a 
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slab of a periodic lattice, and thus, despite being aperiodic, these tilings show a 
pronounced regularity. The numerical results unanimously corroborate the pre- 
diction of the Harris-Luck criterion, and, undoubtedly, these models belong to 
the same universality class as the square-lattice Ising model. 

Apart from these predominantly numerical approaches, also some analyti- 
cal techniques have been employed. This concerns, for instance, high- and low- 
temperature expansions |l9| that can be adapted to quasiperiodic systems. For 
the Penrose and the Ammann-Bccnkcr tiling, the high-temperature expansion 
of the relevant part / of the free energy, 

oo 

/ = E 92n w 2n , w = tanh(/3J) , (30) 

where J denotes the coupling constant, has recently been calculated to the 18th 
order in w JmJ. The coefficients g 2n can be calculated from the frequencies of 
certain graphs of circumference 2n in the tiling, which, for cut-and-project sets, 
can be calculated explicitly. Furthermore, each graph carries a certain weight; 
these weights can be computed recursively. As an example, the expansion for 
the Ammann-Beenker tiling is given by 

w 4 + Xw 6 + (47±-17A)w 8 + (138-50A)ui 10 + (803±-310±A)u; 12 
+(586A-1220)u> 14 + (96§ + 295iA)w 16 + (46566±A-108706)w 18 + 0(w 20 ) 

where A = 1 + y/2. This may be compared with the square-lattice result W9[ 

w A + 2w e + 4±w 8 + 12w w + 37±w 12 + 130w 14 + 490±w 16 + 1958§w 18 -I- 0(w 20 ) 

where the coefficients are rational numbers. From the radius of convergence of 
the series, and the behaviour close to it, one can, in principle, derive the critical 
temperature and the thermal critical exponent. However, the information about 
the critical behaviour that one can extract from the expansion is rather poor, 
because, in contrast to the square-lattice case, the extrapolated values show 
extremely strong fluctuations ]5(],[5l), and much more terms would be necessary 
in order to give reliable estimates of the critical temperature and the critical 
exponents. But this is not feasible; a total of 244 638 different subgraphs of the 
Ammann-Beenker tiling contribute to the coefficient g\%, as compared with a 
mere 1975 different graphs in the square-lattice case. A series analysis of other 
quantities such as the magnetic susceptibility may yield more stringent evidence, 
but this has not yet been calculated because the number of graphs that have to 
be considered is even larger. 

Another method that relies on the computation of the zeros of the partition 
function, in the complex temperature variable z = exp(2J/k B T). The set of 
zeros accumulates on certain curves or areas in the complex plane that separate 
different analytic domains, hence different phases of the system. Thus, a zero on 
the real positive axis will correspond to a phase transition point of the model, 
and the zeros close to this point contain information about the critical exponents. 
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Fig. 5. Part of partition function zeros in the complex variable z — exp(2J/fc B T) for 
a periodic approximant of the Ammann-Beenker tiling with 41 spins per unit cell 

For periodic approximants of aperiodic tilings with rather large unit cells, the 
zero patterns can be calculated explicitly ]5l| . The result for an approximant 
of the Ammann-Beenker tiling, shown in Fig. |5|, is rather involved, whereas for 
the square lattice case, the zeros are restricted to two circles with radius y/2, 
centred at z = ±1. The two intersections of the zero pattern with the real 
axis correspond to the ferromagnetic (J > 0) and antiferromagnetic (J < 0) 
critical point, respectively, which are related to each other because the tiling is 
bipartite. The corresponding estimates of the critical temperature are given in 
Table |l|. Clearly, this method allows a very precise determination of the critical 
temperature, in good agreement with, less accurate, Monte Carlo results ]5^ ]. 

Table 1. Critical temperatures w c — tanh( J/K B T C ) for periodic approximants of the 
Ammann-Beenker tiling, and estimated value for the aperiodic tiling 



m 


N 


Wc 


1 


7 


0.396 850 570 


2 


41 


0.396 003 524 


3 


239 


0.395 985 346 


4 


1393 


0.395 984 811 


5 


8119 


0.395 984 795 


6 


47321 


0.395 984 795 


oo 


oo 


0.395 984 79(1) 
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5.3 Aperiodic Ising Quantum Chains 

We now turn our attention to the case of one-dimensional aperiodicity and con- 
sider Ising quantum chains with aperiodically modulated coupling constants |18| 



H : 





(31) 



where, in contrast to (|l4|), the couplings Ea G {e a ,£f,} an< ^ the transversal fields 
hj £ {h a , h b } are chosen according to the jth letter of an aperiodic substitution 
sequence. 

Such models can be treated by an exact real-space renormalization approach, 
exploiting the recursive structure that is already inherent in the substitution se- 
quence ]53] , ^ , p^y5^ , p7[ . The basic idea behind this approach is that a renormal- 
ization step reverses a substitution step, i.e., the Hamiltonian is transformed 
into a Hamiltonian at the previous substitution level, which just differs from 
the original Hamiltonian by renormalized values of the parameters. The renor- 
malization transformation then becomes a mapping in the finite-dimensional 
parameter space of the Hamiltonian, and the renormalization transformation 
is exact because no additional parameters have to be introduced. This is very 
interesting because it provides one of the few examples with an exact renormal- 
ization transformation for a non-trivial model, keeping in mind that the usual 
decimation approach does not work even for the square-lattice Ising model. Fur- 
thermore, it works for arbitrary substitution sequences |^,^7|, not just for our 
examples (|26| ) , and can even be applied to the case of random substitutions |27| . 
This encompasses a large class of models with different fluctuations in the cou- 
plings and, consequently, different physical behaviour. For the sake of space, we 
cannot go into detail here, but just briefly sketch the calculation. 

The first step consists of a mapping of the Hamiltonian ( |3l| ) onto a model of 
free fermions, 



N 



H = ^2 A j VjVj + C , {r)j,m} = {Vj>Vt} = , {Vj>Vt} = s . 



3,t 



(32) 



by means of a Jordan- Wigner transformation 24 . In this way, the task of com- 
puting the spectrum of (|3|) is reduced from the diagonalization of a 2 N x 2 N 
matrix to that of the 2N x 2N eigenvalue problem 



/ k-y 

h x e-y 
£! 



hn 



h e 9 















••• Oe 
\e n ••• I) // 



#1 



\ -i 



N 



0/ 



A 



N 



= A 



<I>,_ 



* N 



N 



\#N J 



(33) 
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whose N positive solutions < A 1 < A 2 < • ■ • < A N determine the entire 
spectrum of the Hamiltonian ( pl| ) . The quantum phase transition appears when 
the lowest excitation energy tends to zero, A 1 — > 0, in the thermodynamic limit 
N — > oo, so we can concentrate on the low-energy behaviour. 

The renormalization equations can now be derived by exploiting the recursive 
structure of the sequence, effectively combining those sites that belong to a single 
letter at the previous stage of the substitution rule. This results in a mapping 
in parameter space, which can be expanded in powers of A at the fixed point 
A = 0, which corresponds to the critical point given by the condition 



1 = lim 





El 


£ N 


N 




Pa 




K 


h 2 ' 


' h N 




K 




K 



where the last equality applies to the two-letter case discussed above. It can 
be parametrized as \e a /h a \ = r~ Pb and \s b /h b \ — r p - with r S R + . The map- 
ping determines the finite-size scaling of the smallest excitation energy Ai, and 
thus the critical behaviour. In complete accordance with the Harris-Luck cri- 
terion, we find that A\ ~ v{r)N~ x for irrelevant aperiodicity, as in the peri- 
odic case. For marginal aperiodicity, we observe A 1 ~ ]\[- z ( r ) with a coupling- 
dependent non-universal exponent z(r), corresponding to an anisotropic scal- 
ing of the correlation length with different exponents and Uj_. For relevant 
aperiodicity, we find that the lowest excitation energy vanishes exponentially, 
A 1 ~ exp[— Z\(r)7V w ], where ui is the fluctuation exponent (|l9|). In fact, we can 
calculate also the non-universal coupling-dependent terms. For our examples at 
hand, this yields v(r) = 21n(r)/(r — 1/r) for the Fibonacci chain k = 1 (fztf), the 
exponent z(r) — ln(r 1//3 + r -1 / 3 )/ ln(2) for the period-doubling sequence k = 2 

([H, and the coefficient A(r) = v / 2\/l3 - 7\ ln(Y)| for the binary non-Pisot se- 
quence k = 3 J26|). This renormalization approach can also be applied to obtain 
information about the eigenvector in (|33"|), thus determining the behaviour of 
the surface magnetization in an open quantum chain f36| . Furthermore, other 
free-fermion models such as the XY chain can also be treated p7| , |57| | . 



6 Summary and Conclusions 

The article gives a brief overview on recent results for aperiodic Ising models. 
According to the personal taste and scientific experience of the author, aperiod- 
ically ordered Ising models received most attention, while random Ising models 
were only discussed from a qualitative perspective. It is shown that heuristic 
scaling arguments like the Harris-Luck criterion can provide a rather powerful 
tool that, in most cases, correctly predict whether or not a certain disorder will 
affect the critical behaviour. Still, behind the scene are a number of tacit as- 
sumptions, and there exist at least exceptional systems which defy the relevance 
criterion, as shown by an exactly solvable example. 

It is quite a difficult task to treat randomly disordered or aperiodically or- 
dered planar Ising models analytically, and rather sophisticated techniques have 
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been employed. However, it turns out that the Ising quantum chain with cou- 
pling constants modulated according to substitution sequences can be solved by 
an exact real-space renormalization approach, thus proving the validity of the 
Harris-Luck criterion for this entire class of models and, in addition, providing 
information about the non-universal coefficient that enter the scaling behaviour. 
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